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ABSTRACT 

We present the first results of our optical spectroscopy program aimed to provide redshifts and 
identifications for the X-ray sources in the Extended Chandra Deep Field South. A total of 339 
sources were targeted using the IMACS spectrograph at the Magellan telescopes and the VIMOS 
spectrograph at the VLT. We measured redshifts for 186 X-ray sources, including archival data and 
a literature search. We find that the AGN host galaxies have on average redder rest-frame optical 
colors than non-active galaxies, and that they live mostly in the "green valley" . The dependence of the 
fraction of AGN that are obscured on both luminosity and redshift is confirmed at high significance 
and the observed AGN space density is compared with the expectations from existing luminosity 
functions. These AGN show a significant difference in the mid-IR to X-ray flux ratio for obscured 
and unobscured AGN, which can be explained by the effects of dust self- absorption on the former. 
This difference is larger for lower luminosity sources, which is consistent with the dust opening angle 
depending on AGN luminosity. 

Subject headings: galaxies: active ~ quasars: general - X-rays: galaxies 



1. INTRODUCTION 

X-ray surveys have been critical in unveiling a large 
population of obscured s upermassive black ho les actively 
accreting matter (e.g., iGiacconi et al.l [20011 and refer- 
ences therein), contrary to optical su rveys, like th e Sloan 
Digital Sky Survey (e.g., [Schneider et al.ll2002l ). which 
are dominated by unobscured Active Galactic Nuclei 
(AGN). The hard spectral shap e of the X-ray Back- 
ground (XRB; e.g. iGrubeii [19920 reveals that obscured 
AGN should significantly outnum ber the unobscure d 
sources, as was first proposed by Set ti fc Woltieil (|1989D . 
Given the difhculties in detecting and identifying the low 
luminosity and/or high redshift AGN, in particular those 
heavily obscured, it has been very challenging to obtain 
a complete census of the AGN population. This is par- 
ticularly important in order to study the growth of su- 
permassive black holes and its relation to the formation 
and evolution of its host galaxy. 

Deep X-ray surveys like the Chandra Deep Fields 
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North (CDF-N: iBrandt et alllMil and South (CDF-S; 
IGiacconi et al.l 120011 iRosati et al.ll200l have been crit- 
ical in obtaining a significant sample of obscured AGN 
up to significant distances. The first results from the 
optical follow-up of the X-ray sources in the Chandra 
Deep Fields reported a pe ak of the AGN activity at 
(jSzokolv et al.l [2"004[ ) . in contradiction with the 
predictions of the early XRB populat ion synthesis (e.g., 
IComastri eraLlll995l:[Gini et al.ll2001h . which expected a 
redshift peak at This discrepanc y was resolved in 

the XR B populatio n synthesis models of lTreister &: Urrvl 
(|2005[ ) and later bv lGilh et all (|2007D . both of whom in- 
corporated the latest Chandra and XMM-Newton ob- 
servational results in their calculations. Similarly, the 
observed fraction of obscured to unobscured AGN in the 
Chandra Deep Fields was ~1-2:1, lower than the value of 
~3-4:l expected from XRB population synthesis models 
and o bserved in samples of near by sources ( | R ,isalit^ et al.l 
I1999D . This was explained bv iTreister et all ()2004D as 
due to the difficulty in obtaining spectroscopic identifica- 
tions for obscured AGN at significant redshifts, which are 
much fainter in the optical than the unobscured sources. 

One significant problem of these deep multiwavelength 
surveys is caused by the eff ects of cosmic variance. For 
example, iGilli et aLl (|2005f) found a factor of two dif- 
ference between the correlation length for AGN de- 
rived in the CDF-N and CDF-S fields, each covering 
'^O.l deg^, and concluded that this difference could be 
explained by the effects of cosmic variance. In addi- 
tion, significant overdensities and lar ge scale structure s 
have been found in both fields (e.g., iGilh et al]|2003f) . 
which potentially affect the conclusions derived from 
each field individually. In order to address this prob- 
lem, recently several wide-field deep surveys have been 
carried out. Examples of these surveys are the AU- 
waveleng th Extended Grot h strip International Survey 
(AEGIS; iDavis et al ll2007D. t he Cosmic Evolution Sur- 
vey fCOSMOS: IScoville et al]|2007D and t he Multiwave- 
length Survey by Yale-Chile (MUSYC; iGawiser etld] 
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I2006at iTreister et all I2007D . The Extended Chandra 
Deep Field South (ECDF-S) is one of the four 30'x30' 
fields studied by MUSYC. This field is ~3x larger than 
the existing CDF-N or CDF-S. Hence, while the effects 
of cosmic variance are reduced compared to the results 
of each of these deep fields, the uncertainty due to large 
scale structure can still be relevant for our results. To 
minimize these effects, when possible, we combine in 
this work our data with available measurements in other 
fields. 

Deep optical imaging has been obtained in the ECDF- 
S, mostly using the Wide Field Imager (WFI) on the 2.2- 
metre telescope at La Silla. Images in a narrow band fil- 
ter centered at 500 nanometers were also obtained in this 
field in order to look for Lyman Alpha Emitters at z~3 
(jGawiser et al.ll2006bD . The ECDF-S im a.ging was made 
publi c by the ESO D eep Public Surve y (|Mignano et all 
[200l . COMBO- 17 (iWoff et all l200l. and Gar chin 
Bonn Deep Survey (GaBODS: iHildebrandt et al.l[20& 
teams. Deep near-IR coverage was obtained using the 
CTIO 4-metre telescope with the Infrared Sideport Im- 
ager (ISPI), reaching a magnitude limit of ~22 (AB) in 
the JHKs bands^° (E. Taylor et al. in prep.). In ad- 
dition, the ECDF-S is the target of two Spitzer legacy 
surveys, the Spitzer IRAC / MUSYC Pubhc Legacy in 
ECDF-S (SIMPLE) and the Far-Infrared Deep Extra- 
galactic Legacy Survey (FIDEL), which obtained deep 
mid-IR data from 3 to 70 microns. Most relevant for our 
work is the Chandra coverage of the ECDF-S, consisting 
of 4 ACIS-I pointings centered on the original CDF-S to 
a depth of ~230 ksec and covering an area of ~0.3 deg^. 
A detailed de s cripti on of these data was pre sented by 
iLehmer et all (|2005D and IVirani et al] (|2006f ). In this 
work we use the IVirani et al.l ()2006f ) catalog for the op- 
tical follow-up. A total of 651 sources are included in 
this catalog to a flux limit of 1.7xl0~^^ erg cm~^s~^ in 
the 0.5-2 keV band. These sources constitute the main 
sample of this work. 

In this paper we present the first results of an extensive 
spectroscopic effort aimed to identify a significant num- 
ber of the X-ray sources in the ECDF-S field. Optical 
spectroscopy allow us to measure redshifts, and hence 
luminosities, and also to provide an indication of the 
nature of the source, which can be compared with the 
observed multiwavelength properties. In addition, we 
will focus on the properties of the obscuring material 
by studying the mid-IR properties of these sources, and 
by measuring the fraction of obscured AGN and its de- 
pendence on parameters like luminosity or redshift. In 
Section 2 we present the spectroscopic data and out- 
line the reduction techniques used in this work, while 
in §3 the results are presented. In Section 4 we take ad- 
vantage of the exquisite multiwavelength coverage avail- 
able on the ECDF-S to compare the X-ray, optical and 
mid-IR properties of these sources. We also study the 
statistical properties of these sources in order to under- 
stand the evolution and physical conditions of the AGN 
population. The conclusions of our work are presented 
in §5. When required, we assume a ACDM cosmology 
with /io=0.7, r2m=0.3 and r2A=0.7, in agr eement with 
the most recent cosmological observations (Sperg el et al.l 

Optical an d near-IR images and catalogs ar e public and can 
be obtained at |http://www. astro. yale.edu/MUSYC | 



l2007f) . All magnitudes in this paper are pres ented in the 
AB photometric system (|Oke fc Gunnlll983^ . 

2. OPTICAL SPECTROSCOPY 

Sources targeted for spectroscop y were select ed from 
the catalog of Chandra sources of IVirani et all ([2006). 
W e use this catalog in stead of a similar one presented 
bv ILehmer et al.l (|2005( ) since the former uses a more de- 
tailed rejection of high-backg round periods , thus elimi- 
nating spurious sources. The IVirani et al.l (|2006f ) cata- 
log consists of 651 X-ray detected sources, 587 detected 
at high significance (the primary catalog) and 64 found 
using a lower threshold (secondary catalog). The X- 
ray sources were then matched to the MUSYC BVR 
optical catalogue presented by iGawiser et al.l (j2006bl ) 
using a likelihood rat io technique, as described by 
ICardamone et all (|2008f ). using a reliability threshold of 
0.6. The estimated number of false matches using this 
technique, calculated by randomly changing the X-ray 
positions, is ~1%. We found optical counterparts for 445 
X-ray sources in the primary catalog (76%) and 28 in the 
secondary catalog (44%) down to a limiting magnitude of 
BVR=27.1. Using the Spitzer IRAC observations of the 
ECDF-S described in §2.5 and a similar likelihood ratio 
method we found near-IR counterparts for 554 (94%) X- 
ray sources in the primary catalog and 43 (67%) in the 
secondary catalog. Hence, we are very confident that a 
very large fraction of the X-ray sources, >90% in the 
primary and ~70% in the secondary catalogs, are real. 

While in principle all X-ray sources with optical coun- 
terparts were targeted for spectroscopy, we gave higher 
priority to the optically brightest sources, in order to 
increase the efficiency of our observations. In general, 
sources with i? <24 had a higher priority, except for our 
VIMOS (VLT Visible Multi-Object Spectrograph) obser- 
vations, in which an R<25 threshold was used. A total 
of 235 sources with i?<24 and 292 with i?<25 were tar- 
geted. As an experiment, for our IMACS (Inamori Mag- 
ellan Areal Camera and Spectrograph) run of November 
2005 we gave a higher priority to a group of 70 X-ray 
sources with bright optical counterparts (i? <24) and a 
hard X-ray spectrum, given by a hardness ratio (HR) 
lower than -0.2. This was done in order to see if these 
sources have signs of obscuration in the optical spectrum 
as well, and to try to partially overcome the bias against 
obscured sources in identification studies based in opti- 
cal spectroscopy. All these possible selection effects are 
considered on the analysis presented in 521 

2.1. Magellan IMACS data 

We used the IMACS wide-field spectrograph mounted 
on the Magellan I Baade telescope at Las Campanas Ob- 
servatory. This camera provides a circular field of view 
with a 27.2' radius in the //2 configuration. The pixel 
scale for this setup is 0.2"/pixel. For all our observations 
we used the 300 lines/mm grism and slits of 1.2" width, 
except for the October 2003 and 2006 runs in which we 
used 1" slits. For our average seeing of 1", this trans- 
lates into resolution elements of 8 A (6. 5 A for the Octo- 
ber 2003 and 2006 runs) except for the November 2005 
run in which the average 0.6" seeing means a resolution 
element of ~4A for unresolved sources. While the wave- 
length coverage depends strongly on the position of the 
source in the mask, we required a minimum wavelength 
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coveraee of 4200-7000A , which maximizes the number 
of sources on a given mask while at the same time pro- 
viding a broad wavelength coverage, which is important 
considering that the redshifts are not know a priori and 
can span a wide range. The log of our ECDF-S IMACS 
observations is presented in Table [T] 

Typically each mask was observed for 5 hours. The 
goal was to reach a magnitude of R~24, however in many 
cases, due to weather or technical problems, this was not 
possible. Hence, the efficiency, defined as the fraction 
of sources identified and with a measured redshift, dif- 
fers significantly from mask to mask. The total num- 
ber of sources on a typical mask is ~80-100, of which 
15-50% are X-ray sources. The large spread in the frac- 
tion of X-ray sources on a given mask reflects the fact 
that these sources were not always the main target of the 
observations, but "mask-fillers" . The remaining targets 
include extragalactic sources like Ly-a emitters, Lyman 
break galaxies and galactic white and brown dwarf stars 
selected by their optical colors from the MUSYC cata- 
logues. The properties of the general sources targeted 
with IMACS will be presented in a following paper (P. 
Lira et al. in prep.). 

2.2. VLT VIMOS data 



The VIMOS (jLePevre et al.l l2003f ) spectrograph 
mounted at the Nasmyth focus B of UTS was also used 
to obtain optical spectroscopy for the X-ray sources in 
the ECDF-S. This instrument consists of four cameras 
of 7'x8' each, hence providing a total field of view of 
15'xl6', with a pixel scale of 0.2"/pixel. We used the 
MR grism together with the OS-blue filter in order to 
block the contamination from higher orders. In this setup 
the wavelength coverage is 4000-6700A , very similar to 
the coverage of the IMACS observations. This coverage 
minimizes the effects of fringing at red wavelengths and 
takes advantage of the wavelength region of maximum ef- 
ficiency of the VIMOS camera. The MR grism provides 
a dispersion of 2.5A /pixel, which combined with our 1" 
wide slits corresponds to a resolution element of 12. 5A . 
All the VIMOS observations were carried out in service 
mode by the Paranal staff. The required conditions were 
seeing better than 1", clear sky and fractional lunar illu- 
mination smaller than 20% at a minimum distance of 30 
degrees. Due to the lack of an Atmospheric Dispersion 
Corrector in VIMOS and in order to minimize the effects 
of instrument flexures, observations were carried out at 
a maximum hour angle of ±2 hours. 

Four pointings were used to completely cover the 
ECDF-S 30' X 30' fleld. Each mask was observed for three 
hours and 18 minutes, which corresponds to a total exe- 
cution time of six hours including overheads. According 
to the exposure time calculator this should allow us to 
significantly detect sources up to R^24.5. A total of 283 
sources were targeted, excluding alignment stars, 96 of 
them X-ray sources. 

2.3. Additional Observations 

Since the ECDF-S is the target of many deep multi- 
wavelength observations, it is not surprising that many 
spectroscopic follow-up programs have been carried out 
in this field. In order to take advantage of existing ob- 
servations of these X-ray sources, we did an extensive 



literature search. In particular, we used the Master Com- 
pilation of GOODS/CDF-S spectroscopy^^. We matched 
this spectroscopic catalog to our MUSYC optical cata- 
log using a search radius of 0.7". Then, we used the 
MUSYC-Chandra matched catalog described before and 
found spectroscopic identification s for the X-ray s ources 
in the obs e rvations reported by IVanzella et al.l (I2005L 
[2006'.'2008').'Le Fevre et"all ((200l T lSzokolv et al.l (l200l ) 
and Croom et al. (2001). In total, existing identifications 
of 55 X-ray sources were added to our sample. 

We also searched for existing unpublished spectro- 
scopic observations in the ESO-VLT archive^^. We fo- 
cused on VIMOS, given the large area covered by this 
camera and that most of the existing F0RS2 observa- 
tions were already incorporated in the master compi- 
lation described above. There are several VIMOS ob- 
servations of the ECDF-S on the ESO archive, however 
of particular relevance for us was the program of Berg- 
eron et al. (program ID 072.A-0139), for which 23 hours 
were granted in order to obtain spectroscopy of the X-ray 
sources detected by XMM. The area covered by XMM is 
well matched to the ECDF-S Chandra observations, so 
it is reasonable that a significant number of our Chandra 
sources were identified as part of this program. Three 
fields covered by the Bergeron program overlap with the 
ECDF-S. These fields were observed using both the MR 
and LR-Blue grisms with the GG475 and OS-Blue filters 
respectively. We looked for counterparts of our Chandra 
sources using a search radius of 2", which was chosen in 
order to account for the XMM positional uncertainty but 
at the same time avoid mis-identifications. In total, 36 
sources in our Chandra catalog were observed as part of 
the Bergeron program; 28 of them were already targeted 
as part of our IMACS and VIMOS program, but only 11 
successfully identified, so potentially identifications for 
25 new sources could be obtained from the archival data. 



2.4. Data Reduction 

The IMACS data were reduced using the Carnegie Ob- 
servatories System for MultiObject Spectroscopy (COS- 
MOS) package^^ developed by Gus Oemler. This soft- 
ware was specifically designed to reduce IMACS data, 
in particular handling the challenge of having spectra 
spread over several CCDs. We followed the standard re- 
duction cookbook, which consists of the following steps: 
construction of the spectral map using arc frames, bias- 
subtraction and flat fielding of science data and sky 
subtraction. We used COSMOS to generate 2D sky- 
subtracted wavelength-calibrated spectra of each source. 
We then used IRAF's apall task in order to extract the 
ID spectrum. 

To reduce the VIMOS data we used the ESO VIMOS 
pipeline version 2.1.6^'*. In this case all the spectra for 
a given quadrant are located on one CCD, so the data 
reduction is significantly simpler than for IMACS data. 
However, the basic reduction steps are very similar to 
the IMACS reduction. The same procedure was used 

Available at http:/ /www. eso.org/science/goods/spectroscopy/CDFS_Maste 
The ESO archive can be found at http://archive.eso.org] 

This software can be obtained from 



I http: //www.ociw^e du / Code /cosmos| 

This package can be obtained from 

|http://www.eso.org/projects/dfs/dfs-shared/web/vlt/vlt-instrument-pipelines.h 
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to reduce both the archival observations and our pro- 
prietary data. In the case of the VIMOS data from our 
program, for each individual observation of ~33 minutes, 
the corresponding night-time arc was used to calculate 
the spectral map, in order to account for the effects of 
instrument flexures. In addition, observations of stan- 
dard stars on the same night were used to calculate the 
effective response as a function of wavelength for these 
observations. 

2.5. Spitzer Data 

In order to study the IR to X-ray ratio for our sam- 
ple of sources with measured spectroscopic redshifts in 
the ECDF-S we took advantage of the deep Spitzer im- 
ages available in this fleld. At the shortest wavelengths, 
the ce ntral CDF-S regi on (10'xl6') was observed by 
IRAC (jFazio et al.ll2004D as part of the GOODS survey 
to flux limits of 0.13, 0.22, 1.45 and 1.61 ^Jy in the 3.6, 
4.5, 5.7 and 8 /xm band respectively. So urce matching 
and b asic IR properties were reported by iTreister et al.l 
(Hooe"). In addition, the whole ECDF-S field was cov- 
ered by IRAC as part of the SIMPLE survey^^ (Damen 
et al. in prep.). The flux limits for the SIMPLE ob- 
servations are 0.76, 0.4, 5.8 and 3.6 /iJy in the 3.6, 4.5 
and 5.7 and 8 ^m bands respectively, thus '^3-5 x shal- 
lower than the GOODS observations. The matched cat- 
alog of the ECDF-S X-ray sources to the SIMPLE im- 
ages used here was pres ented and described in detail by 
ICardamone et al.l 



At l onger wavelengths , we use the available Spitzer 
MIPS (jRieke et al.| l2004^ observations at 24 /im. In the 
central region, GOODS provide deep images to a flux 
limit of ~25 /iJy. A c atalog^^ based on PS F-fitting mag- 
nitudes, described bv lTreister et al.l (|2006[ ) and R. Chary 
et al. (in prep.), was matched to the X-ray sources in the 
GOODS region using the IRAC-derived positions. 

Additional data at 24 /im covering the whole ECDF- 
S field were taken as part of the FIDEL survey. Un- 
fortunately, at the moment only the reduced images 
were made available by the FIDEL group^^. In order 
to produce our own catalog based on the reduced FI- 
DEL images we used the Mosaicking and Point- Source 
Extraction fMOPEX: fMakovoz fc Marleaul [2005[) pack- 
age, provided by the Spitzer Science Center. Briefiy, 
we did a flrst-pass extraction using the default point 
response function (PRE) included in MOPEX in order 
to obtain a preliminary catalog of the brightest sources. 
These sources were then used to construct a customized 
PRF, based directly in the FIDEL observations. Fi- 
nally, using this PRF we performed a second extraction 
and measured PRF-fitting flux densities for the detected 
sources. This procedure is described in more detail by 
iMakovoz fc Marleaul ()2005( ). In order to verify the ac- 
curacy of our derived fluxes in Figure [T] we compare 
the FIDEL and GOODS sources for the 829 overlapping 

Further information about this survey, cat- 
alogs and images can be f ound at the website 
|http: / /www. astro. yale.edu/dokkum/simple/ 

^° More details about Fhis catalog and data- 
reduction procedures can be found at the webpage 




[http: / / data. spitzer. caltech.edu/popular / goods/Documents / goods. 

^' The reduced images and more details about 
the survey and reduction techniques can be found at 
|http://data.spitzer. caltech.edu/popular/fidel/2007-sepl7/fidel-dr 2 



GOODS flux (/i,Jy) 

Fig. 1. — Ratio of PSF-fitting 24 /^m fluxes derived from the 
FIDEL images and GOODS catalogs as a function of GOODS fiux 
for the sources significantly detected by both surveys. The median 
value for the flux ratio is 0.97, with a standard deviation of 7%, 
well within the measurement errors, showing that the two studies 
produced consistent results. 

sources with fluxes higher than 60 /iJy in order to en- 
sure a reliable detection in the FIDEL image. There is a 
very good agreement between the FIDEL and GOODS 
fluxes. The median ratio of FIDEL to GOODS flux is 
0.97 with a standard deviation of 7%. Given the very 
small systematic offset between the FIDEL and GOODS 
fluxes, well within the measurement errors, no correction 
to the FIDEL fluxes was applied. The FIDEL catalog 
was matched to the ECDF-S using the maximum likeli- 
hood procedure described by Carda mone et al.i f)2008) . 

3. RESULTS 

A total of 339 X-ray sources in the lVirani et all (|2006f ) 
catalog were targeted for spectroscopy using the VIMOS 
and IMACS spectrographs and data from the literature. 
For 186 of them, we were able to identify the nature of 
the source and measure the redshift, giving an average 
efficiency of 55%. Of those, 180 were detected in the 
high-signiflcance primary X-ray catalog. The basic X- 
ray, optical and spectroscopic properties of these sources 
are presented in Table [2] The R-band magnitude dis- 
tributions for all the X-ray sources, the sources targeted 
for spectroscopy and those with spectroscopic identifica- 
tion are presented in Figure [2] The fraction of sources 
that are successfully identified depends strongly on the 
observing conditions, achieved exposure time, etc. (Ta- 
ble [T]). At the same time, this fraction always decreases 
strongly towards fainter magnitudes (Figure [2]) . In our 
analysis, we model this dependence using a very simple 
linear fit, with a constant value of 65% to R=22 and then 
decreasing linearly to 0% at i?=26.5. The resulting fit is 
shown in the upper panel of Figure [2l 

The redshift distribution for our spectroscopically- 
identified sources is presented in Figure [31 The vast ma- 
ioritj^ of unobscured sources (that present broad emission 

datap^o''"-**" I-*™' 
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fines, see discussion below) are at z>l, with an average 
redshift of 2.22 and a median of 2.12. In contrast, ob- 
html scured AGN and X-ray emitting galaxies are preferen- 
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Fig. 2. — Bottom panel: Distribution of ij-band magnitude for 
all X-ray sources with optical counterpart in the MUSYC images 
{solid line), sources targeted for spectroscopy (dashed) and identi- 
fied sources (hatched histogram). Top panel: Fraction of identified 
to total X-ray sources as a function of R-band magnitude. The thin 
solid line shows the fit to the data points, used in §4. The frac- 
tion of identified sources is ~65% at R<22, and decreases towards 
fainter optical fiuxes. 




Redshift 

Fig. 3. — Bottom panel: Redshift distribution for all X-ray 
sources with spectroscopic identification (solid line) and sources 
classified as unobscured AGN only (hatched histogram) in bins of 
A2=0.2. Top panel: Redshift distribution for all sources in bins 
of Az=0.02, in order to identify overdensities. The dashed lines 
mark the previously known structures in this field at 2=0.67 and 
0.73, while the dotted lines show other suspected overdensities at 
2=0.15, 0.29 and 0.53. 

tially located at z<l. For the total sample, the average 
redshift is 1.18, with a median of 0.75, in good agree- 
ment with the result s found in both Chandra Deep Field s 
(|Szokolv et al.l[200l iBareer et al.l[200l iGilli et aLllMsl) . 
The highest redshift source in our sample is a broad line 
AGN, XID 213, at z=4.48. 

In the upper panel of Figure [3] we use a narrower bin 
width (Az=0.02), in order to search for structures in red- 



FlG. 4. — Solid small circles: location of identified sources in 
the ECDF-S field. Other symbols mark the sources in the identi- 
fied overdensities: Solid triangles: z=0.15, solid sguares: 2=0.29, 
large solid circles: 2=0.53, open squares: 2=0.67, open triangles: 
2=0.73. As can be seen, sources in the groups in redshift space 
span almost the whole studied area. 

shift space. The two main structures, "walls" , previously 
detected in. the CDF-S proper at z=0.67 and 2:=0.73 
(|Gim et al.ll200llAdami etal][200l are also easily found 
in our sample. In addition, we found excesses of sources 
at z=0. 15,0.29 and 0.53. In Figure|4]we present the spa- 
tial distribution of all the identified sources, with the 
position of the sources in those redshift structures high- 
lighted. For z=0. 15,0.29 and 0.53 the 30'x30' ECDF-S 
field corresponds to a comoving area of 4.6x4.6, 7.7x7.7 
and 11.2x11.2 Mpc^ respectively. Given that all the red- 
shift groups span almost the whole field of view, these 
structures a re much large r than the compact structures 
reported bv lAdami et al.l (|2005f ). In order to estimate 
the probability of having seven sources in one Az=0.02 
bin we assume that the redshift distribution follows a 
Poisson distribution. The average number of sources in 
a Az=0.02 bin is ~2. Hence, the probability of hav- 
ing seven sources in one bin randomly is ~0.3%. There- 
fore, we can conclude that these groups in redshift space 
should correspond to real physical associations. 

Sources were mainly classified according to their op- 
tical spectroscopy characteristics. Our primary division 
comes from the presence or absence of broad emission 
lines. The threshold for unobscured sources was set at 
Av>1000 km s~^, since Bargcr et al. (2005) showed that 
the vast majority of the non-X-ray sources have line- 
widths smaller than this value, while sources with broad 
lines have also soft X-ray spectra, consistent with be- 
ing dominated by an unobscured AGN. We then used 
the hard (2-8 keV) X-ray luminosity to separate nor- 
mal galaxies with X-ray emission from AGN-dominatcd 
sources. Specifically, we assumed a conservative and typ- 
ical threshold of Lx=10^^erg s^^, higher than the high- 
est level of X-ray activity observed in star-forming galax- 
ies (e.g. JLira et an[200l . 

The X-ray spectral shape can be used to obtain an 
alternative source classification. In particular, the hard- 
ness ratio, defined as the ratio between the difference of 
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Fig. 5. — Hardness ratio, defined as (H-S)/(H+S) where H and S 
are the hard (2-8 keV) and soft (0.5-2 keV) X-ray count rates, as a 
function of hard X-ray luminosity for sources optically classified as 
obscured (open circles) and unobscured [filled circles) AGN. Ver- 
tical dashed lines show the typical separation for normal galaxies 
{Lx<W^^ erg s~^), Seyfert galaxies {Lx=W^^~'^'^ erg s~^) and 
quasars (Lx>10^'* erg s~^). The dashed horizontal line shows the 
HR for a source with r=1.9 and Af^f =10^^cm~2 at z<0.5. Most 
sources classified as unobscured AGN are clustered in the high lu- 
minosity and small HR region, while obscured AGN have in general 
lower luminosities and higher HR. 

the hard band (typicaUy 2-8 keV for Chandra observa- 
tions) and soft band (0.5-2 keV) count rates and their 
sum is used when the number of detected counts is not 
enough to perform X-ray spectral fitting. Figure [5] shows 
the hardness ratio as a function of luminosity for the 
186 sources with spectroscopic identifications presented 
in this paper. As can be seen, sources with hard X-ray 
spectrum (positive hardness ratio) tend to be also opti- 
cally classified as obscured sources in the optical, while 
sources that present broad emission lines have in general 
a soft X-ray spectrum (negative hardness ratio) . It is also 
interesting that most of the unobscured sources, in both 
X-rays and optical wavebands, have a high X-ray lumi- 
nosity {Lx >10*^ erg s~^), while most obscured sources 
are fainter in X-rays. These trends will be studied in 
further detail in the following section. 

While the hardness ratio can be used as a rough in- 
dication of the X-ray spectral shape, for the brighter 
X-ray sources it is possible to perform more accurate 
spectral fitting. In the case of the ECDF-S, a total of 
184 sources have more than 80 counts detected in the 
0.5-8 keV Chandra band, making spectral fitting possi- 
ble. We used a modified version of the Yaxx^^ software 
to extract the spectrum of each high-significance source. 
For sources with more than 200 background-subtracted 
counts and measured redshifts, 72 in the ECDF-S, we 
simultaneously fitted a power-law spectrum and photo- 
electric absorption, with three free parameters: slope 
of the power law (F), normalization and observed neu- 
tral hydrogen column density {Nh)- In all cases we 
included the Galactic absorption, assuming an average 
Nh of 6.8 xlO^^ cm^^ for the ECDF-S, as measured by 

Available at http:/ /cxc. harvard. edu/contrib/ya^cx/] 



the survey of Galactic HI of iKalberla eTall (|2005l ). For 
sources with less than 200 but more than 80 counts it 
was not possible to perform such detailed fitting, so we 
fixed the slope of the power-law to a value of r=1.9, very 
similar to the average value for the sources with a fitted 
slope (r=1.95), which also correspon ds to the typical X- 
ray spectrum for u nobscured AGN (|Nandra fc Poundl 
11994 iNandra et al.|[T997l and others). In both cases, the 
Nh value was derived in the observed frame, in order to 
keep the X-ray spectral fitting independent of the red- 
shift measurements. In order to calculate the conver- 
sion from observed to i ntrinsic Nh w e used the Xspec 
version 12.4.0 software (|Arnaudlll996l ) to simulate a X- 
ray spectrum with r=1.9 as observed at the Chandra 
aimpoint with ACIS-I, added absorption at varying red- 
shifts from to 4 and measured the fitted Nh value 
in the observed frame. We obtained an empirical rela- 
tion of NH{intr)={l+zy -^^NH{ohs), consistent with the 
conversion assume d in the past for Chandra data (e.g., 
iBauer et al.ll200^ . Furthermore, we checked that the 
derived relation is roughly independent of the assumed 
slope of the power-law spectrum. Hence, we are very 
confident that our conclusions are not affected by our 
choice of measuring absorption in the observed frame. 

The conversion from observed to intrinsic Nh can thus 
only be done for sources with a measured redshift. We 
have spectroscopic redshifts, given in Table O for a to- 
tal of 80 sources with fitted Nh- However, in order to 
increase our sample we used the good quality photomet- 
ric redshifts derived by the COMBO- 17 survey, which 
provides red shifts accurate to 5z/{l+z)^10% for sources 
with i?<24 (jWolf et al.ll2004[ ). This uncertainty in red- 
shift translates into an uncertainty in the derived intrin- 
sic Nh of 6\ogNH—0.12, smaller than the typical er- 
rors in the fitted observed Nh or the bin size in our 
Nh distribution (Fig. [6|). Hence, it is acceptable to use 
the COMBO- 17 photometric redshifts in order to con- 
vert observed Nh into intrinsic values, thus increasing 
our sample to 144 sources. To confirm that the use of 
COMBO-17 redshifts does not bias our conclusions, we 
performed a Kolmogorov-Smirnov (KS) test and found 
that the hypothesis that both distributions were drawn 
from the same parent distribution can only be dismissed 
at the 0.02% confidence level. 

4. DISCUSSION 

4.1. X-ray/Optical Classification 

The bottom panel of Figure [6] shows the distribution of 
intrinsic Nh separately for sources classified as obscured 
and unobscured AGN based only on optical spectroscopy. 
While the average value of Nh is larger for obscured than 
unobscured sources, 8x10^^ cm~^ versus 2x10^^ cm~^, 
the difference is not very large. Performing a KS test 
to these two distributions rules out the hypothesis that 
they were drawn from the same parent distribution with 
a 98.5% confidence, leaving a small but not completely 
negligible probability. Since they are both indications 
of the X-ray spectral shape, there is a close relationship 
between the hardness ratio and derived Nh values. As 
can be seen in Figure [71 for modest column densities, 
A^ff<10^^'^ the Nh values are very sensitive to small 
changes in the number of detected soft X-ray photons, 
in particular for high redshift sources. This implies that 
there can be large errors in the derived Nh values, in 
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Fig. 6. — Top panel: Distribution of for all the sources 
with either spectroscopic or photometric redshift and more than 
80 counts in X-rays. The hatched histogram shows the distribu- 
tion for sources with spectroscopic redshifts only. According to a 
KS test, the distributions are entirely consistent with being drawn 
from the same parent distribution. Bottom panel: Nh distribu- 
tion for sources optically classified as obscured {solid histogram) 
and unobscured (dashed histogram). While X-ray obscuration is 
present in unobscured AGN, most of them have very small Nh 
values. Similarly, sources classified as obscured AGN in the op- 
tical have in general higher Nu values. However, in both cases 
significant scatter is present. 
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Fig. 7. — Neutral hydrogen column density, Nh, as a function of 
hardness ratio. Sources optically classified as unobscured AGN are 
shown by solid circles, while obscured AGN are marked by open 
circles. Sources with more than 200 X-ray counts detected are en- 
closed in squares. The dashed lines show the relation between Nh 
and HR for a power-law source with r=1.9 at z=0,l,2 and 3. At 
Af/f <10^^'^, and in particular at high redshift, the Nh measure- 
ments are very sensitive to small changes in HR. This can explain 
the presence of optically unobscured sources with high Nh values. 

particular for unobscured sources at high redshift. The 
imphcations of this potential bias in our analysis will be 
further discussed below. 
As was shown in Figures [5] and [7l the optical and X- 





FlG. 8. — Left panel: Nh as a function of hard X-ray luminos- 
ity for sources optically classified as unobscured (light gray circles) 
and obscured (black) AGN. The large circles show the AGN in our 
sample, while the small circl es are sources in the deeper central 
region with A'^^f measured bv lTozzi et al.l I I2006I) . The lower panel 
shows the luminosity distribution for the same sources. Large sym- 
bols with errorbars show the average values in bins of 20 sources. 
A clear separation in the Nh average can be seen for obscured and 
unobscured sources, except at low luminosities, where the effects 
of dilution by the host galaxy can hide the broad emission lines 
used to classify sources optically. Right panel: Same as left panel, 
but showing Nh as a function of redshift. In this case there is a 
clear increasing trend in the average values with redshift. This can 
be explained by the redshift correction in converting observed into 
intrinsic A'^^f values. Because of this correction, sources optically 
classified as unobscured at high redshift, can appear to have large 
Nh values, as discussed in the text. 

ray classifications do not always agree. In order to fur- 
ther investigate this effect, in the left panel of Fig. [8] 
we show the intrinsic Nh versus hard X-ray luminosity, 
used as a tracer of the intrinsic AGN luminosity. For this 
figure we also added 335 X-ray sources detected in the 
1 Msec. CDF-S observations with accurate Nh measure- 
ments based on spectral fitting presented bv lTozzi et al.l 
(|2006f) . Redshifts for those sources , either phot ometric 
or spectroscopic, were compiled bv IZheng et al.l |2004). 
This figure shows two interesting effects. At low lumi- 
nosities, Lx<W*^ erg s~^, a large fraction of the sources 
optically classified as obscured AGN have intrinsic Nh 
values lower than 10^^ cm~^. This can be explained by 
the ef fects of dilution of the AGN light by the hos t galaxy 
(e.g.. iMoran et all 1200^ ICardamorie et al.ll2007l ). which 
lowers the equivalent width of the broad emission lines 
because of the relatively high continuum from the host 
galaxy, thus making a source to be erroneusly classified 
as obscured AGN in the optical. A second effect ap- 
pears at high luminosity: sources optically classified as 
unobscured AGN have intrinsic Nh values higher than 
10^^ cm~^. A possible explanation for this effect can 
be found in the right panel of figure El where we can 
see that the sources with discrepant X-ray /optical clas- 
sifications are predominantly located at z>2.5. Given 
the strong redshift dependence of the accuracy of Nh 
determinations, sources with very small observed Nh 
values, which can be in general explained by measure- 
ment errors, will have a high intrinsic Nh values if the 
source is located at high redshift. At high redshift, the 
Chandra observations trace emission at higher energy in 
the rest frame, thus making it much harder to detect 
the photoelectric absorption cutoff, critical for measur- 
in g Nh- This e f fect w as already studied and simulated 



bv lAkvlas et al.l (|2006D . who found that at high redshift 
the derived column densities are sy stematically oyerest i- 
mated. According to the results of lAkvlas et al.l ([2006) , 
Nh measurements are systematically overestimated by 
-50% at z-2.5 and -20'% at z-1.5. 
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Fig. 9. — Same as Fig. E] but using the combined classification 
scheme described in the text. In this case, it can be seen that the 
discrepancy between the optical and X-ray classifications for low 
luminosity sources is solved and our hybrid classification clearly 
separate sources in their average Njj values. 

An obvious conclusion of this analysis is that no classi- 
fication method is perfect and can be used for all sources. 
An early attempt to use a combined X-ray /opt ical classi- 
fication was performed bv lSzokolv et al,l ()2004[ ), In their 
work, they used a separation between obscured and un- 
obscured sources at a hardness ratio o f -0.2, A simi- 
lar approach was used by Ifasi nger et al.l (|2005l ) to select 
unobscured AGN from an X-ray selected samples. This 
threshold is shown in Figure [5] for the sources with spec- 
troscopic identification in the ECDF-S. A hardness ratio 
of -0.2 corresponds to an intrinsic Nh of 10^^ cm^^ at 
z=0.5, assuming r=1.9, while sources with this Nh at 
lower redshift will have a higher (more positive) value. 
So, a threshold of -0.2 in HR to separate obscured and 
unobscured AGN is appropriate at z<0.5. However, 
as is clear from Figure [SI the effects of dilution by 
the host galaxy are important only for low luminosity 
sources, which are predominantly located at low redshift. 
Hence, we propose the following classification scheme: 
For sources at z<0.5 we will separate obscured and un- 
obscured sources at HR=-0.2, regardless of their optical 
properties, while for sources at higher redshifts we will 
use only the optical classification scheme, as Nh mea- 
surements are biased at high redshift. 

In Figure [9] we show the measured Nh values as a 
function of luminosity and redshift for the new classifica- 
tion scheme. As expected, changes are relevant only for 
sources with low luminosities at low redshifts. With an 
optical classification only, the average Nh for obscured 
sources at z<0.5 is 4,8x10^^ cm~^, while using the X- 
ray classification it is 7,3x10^"'^ cm~^. Similarly, in the 
lowest luminosity bin, Lx<3,2xl0''^ erg s~^, obscured 
and unobscured AGN are now clearly separated in Nh- 
This separation is not even larger since heavily-obscured 
low-luminosity sources are preferentially missed by the 
X-ray observations, and thus not included in our sample, 

4,2, Rest-frame Optical Colors 

As is well known (eg,. Bal drv et ani2004l : rWeiner et al,l 
l2005t Cirasuolo et al, 2007 and references therein), the 
distribution of optical colors in normal galaxies is bi- 
modal, with a large population of blue star-forming 
galaxies in the blue cloud separated from the "red se- 
quence" of massive pas sively-evo lving spheroids. As was 
recen t ly reported ('e,g,. iBarger et al, 2003; Nandr a et al,l 
I2007t 'Georgakakis et all l2008t ~ ^ilverman et all l2008f ). 
most galaxies hosting obscured AGN, in which the inte- 
grated optical light is dominated by the host galaxy (e,g,, 



iTreister et al,l [20051 ) . lie on the "green valley", a transi- 
tion region located between the red sequence and the blue 
cloud. Hence, this can be considered as good evidence 
for a link between AGN and galaxy evolution, possibly 
indicating that A GN feedback can play a role in regulat- 
ing st ar formation (jSpringel et al.l[2005tlSchawinski et al.l 
120061): however, this reinains c ontroversial, as discussed 
bv lAlonso-Herrero et all (|2008[) . 

To compute the rest-frame optical colors for the X-ray 
sources with spectroscopic identification on the ECDF- 
S field we sta r ted fro m the MUSYC optical catalog of 
iGawiser et al.l ()2006b( ). In order to increase the num- 
ber of sources in our sample we added the X-ray sources 
with reliable photometric redshifts from the COMBO- 
17 survey, as described before. We then generated the 
rest-frame spectral energy distribution for each X-ray 
source with measured redshift by interpolating the ob- 
served UBVRIzJKs photometric data points. The rest- 
frame U and V optical magnitudes were then computed 
by convolving the energy distribution with the corre- 
sponding filter response. The resulting rest-frame U-V 
colors versus My, the rest-frame V^-band absolute mag- 
nitude, are shown in Fig. IIOI A significant number of 
X-ray sources have U-V<1 and My<-22, In general, 
those sources were classified as unobscured AGN because 
of the presence of broad emission lines. Since the opti- 
cal continuum of these sources is dominated by the AGN 
and not the host galaxy, we exclude t hem from further 
analys is. However, as it was found by ISchawinski et aP 
once the emission from the central source is re- 
moved, the host galaxies of unobscured AGN have opti- 
cal colors similar to those of unobscured AGN, As can be 
seen by comparing with a population of normal galaxies 
obtained on the same field from the COMBO-17 survey, 
the X-ray sources classified as obscured AGN have in 
general redder optical colors than normal galaxies and 
lie outside of the main blue cloud. 

The excess of red sources in the AGN population can 
be easily seen in Fig, [Til where the distributions of rest- 
frame U-V colors for X-ray emitting sources and normal 
galaxies are presented. In order to limit the influence of 
differences in the sample selection, the comparison sam- 
ple was obtained randomly from the EGDF-S inactive 
galaxies to match the AGN redshift distribution. Per- 
forming a KS test we found that the probability that 
these two distributions were drawn from the same par- 
ent distribution is negligible, Contrarily, a KS test com- 
paring the distribution of the X-ray sources with spec- 
troscopic and photometric redshift returned a relatively 
high probability of ~20%, indicating that the use of pho- 
tometric redshifts does not significantly bias our results. 
The fact that results of the KS test are not even higher 
can be explained since in general sources with spectro- 
scopic redshifts have brighter optical counterpart, and 
thus slightly different average U-V colors. However, this 
is a minor effect that will not affect our conclusions. 

Given that the comparison sample was selected to 
have the same absolute optical magnitude and redshift 
distributions as the AGN host galaxies, it is unlikely 
that the differences in optical colors are due to the 
sample selection, A sim ilar analysis was performed by 
iGeorgakakis et al.l (|2008 ^ using the AGN detected in 
the AEGIS survey. Studying the morphology of AGN 
host galaxies, they found only a small fraction of major 
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Fig. 10. — Rest frame U-V color versus rest-frame V-band ab- 
solute magnitude for obscured {red circles) and unobscured AGN 
{blue circles). The green boxes mark the location of sources with 
spectroscopic identifications. Normal galaxies (i.e., non-AGN) with 
spectroscopic redshifts are shown by black circles, while the gray 
background shows the density of normal galaxies with photometric 
redshifts from COMBO-17. Contours enclose regions with densi- 
ties larger than 10, 20 and 30 galaxies per interval. While the 
unobscured AGN are in general bluer and have higher luminosities 
than normal galaxies, obscured AGN have similar optical luminosi- 
ties but have redder U-V colors and live arou nd the "gre e n vall ey" , 
shown by the dashed line, as determined by IBell et aP ((2003) for 
z=l. 
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Fig. 11. — Normalized distribution of U-V values for normal 
galaxies {dashed line), AGN {black solid line) and AGN with spec- 
troscopic identification {gray line) for sources with -23<Mv<-18. 
The sample of normal galaxies was selected randomly to match the 
AGN redshift distribution. The separation between normal galax- 
ies and AGN in their U-V colors can be seen clearly. In this range 
of optical luminosities, most of the AGN are obscured, hence the in- 
tegrated emission is dominated by the galaxy, thus indicating that 
in general AGN host galaxies have redder color than non-active 
galaxies. 



mergers, but evidence of minor interactions in a large 
number of source s . Bot h iGeorgakakis et al] ()2008D and 
iSchawinski et al.l (|2008[ ) concluded that AGN activity 
should outlive the truncation of star formation, in or- 
der to explain the observed red optical colors of AGN 
host galaxies. 

4.3. Obscured AGN Fraction 

The fraction of obscured AGN and its possible de- 
pendence on parameters like redshift or luminosity is 
very important for AGN population studies, in partic- 
ular given its direct relation with our understanding of 
the cosmic X-ray Background. A dependence of the ob- 
scured fraction on lumi nosity was previous l y found for X- 
ray se l ected samples bvlUeda et a l. (2003), Stcffcn et alj 
(|2003f l. lLa Franca et al.l (|2005l ) and Treister et ah (,2005f ) 
among others. A physical exp lanation, propo sed by 
iLawrenci (|1991[ ) and updated by ISimpsonI (I2005D. is the 
so-cal led "receding torus" . More recentlv. iTreister et al.l 
([2008') found a dependence of the ratio of mid-IR to 
bolometric flux on luminosity for unobscured AGN con- 
sistent with an incre ase in the opening angle w i th lu - 
minosity. Similarly, lAkvlas fc Georgantopouloi (|2008l ) 
proposed that the observed dependence of the obscured 
fraction on luminosity could be explained by the effects 
of photo-ionization on t he X-ray obscuring matter, while 
IHonig fc BeckCTtI (120071) argues that the effects of the Ed- 
dington limit on a clumpy torus can explain the observed 
luminosity dependence of the obscured AGN fraction. 
The large sample of X-ray sources in the ECDF-S can 
be used to further explore this luminosity dependence 
of the obscured AGN fraction. In figure [12] we can see 
that using the ECDF-S sample alone the fraction of ob- 
scured AGN decreases from ~90% at Lx<10*'^ erg s~^ 
to '--^20% at Lx—W*^ erg s~^, using the "hybrid" classi- 
fication scheme described in §4.1. 

In order to increase the significance of this result, we 
added the res ults from the ECDF-S to the compilation 
presented by ITreister &: Urrvl (|2006f ). which combined 
the results from seven X-ray surveys ranging from wide 
area shallow surveys to the Chandra deep fields. The to- 
tal sample includes now 2814 X-ray sources, 1377 (49%) 
of them with spectroscopic identification. In Figure [iSj 
we present the resulting obscured AGN fraction as a 
function of luminosity for the total sample. The ob- 
scured fraction ranges from 80±7% at Lx<W^^ erg s""'^ 



to 16±8% for Lx>10 



45 



erg s~ . In addition. Fig- 
ure [12] shows the predicted luminosity d ependence in- 
corporated in the ITreister fc Urrvl (I2005D AGN popu- 
lation synthesis model, as updated bv ITreister fc Urrvl 
((2006) to include a redshift dependence, and for com- 
pa rison the luminos ity dependence used in the models 
of iGilh et all (|2007D . While at low luminosities both 
models agree well with the observati ons, for luminosi - 
ties higher than Lx^W^"^ erg s"! the iGilh et "all ()2007f l 
models predict a fraction of obscured AGN of ^50%, 
significantly higher than the observed value. This dis- 
crepancy has important consequences for the modeling of 
the X-ray background using AGN, as the largest contri- 
bution comes from sources at roughly these luminosities 
Treister & Urry 2005) . The larger fraction assumed 



bv lGilli et al.l (|2007i, ) is more relevant at lower energies, 
where the effects of absorption are lar ger. For example , 
the XRB intensity at E<10 keV in the lGilli et al] (|2007f ) 
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Fig. 12. — Ratio of obscured to total AGN as a function of hard 
X-ray luminosity. The black circles with dashed error bars show the 
obscured AGN fraction from the ECDF-S alone, while the gray cir- 
cles show the results obtained using the sample of Treister & Urrv 
12006). Black circles with solid error bars show the fraction of 
obscured AGN for the combined sample of Treister & Urrv ( 200(1) 
and th e ECDF-S sources. The dependence used by the Gi lli et aLl 
II2007I) AGN population synthesis model for the intrinsic and ob- 
served fractions of obscured AGN are shown by the dashed and 
dotted lines. T he dependence used in the model presented by 
ITreister fc Urrvl | |200Q1 is shown by the solid line. While a good 
agreement i s found betwee n the observations and the expectations 
of Treister Sz Urrvl ll2006h . the IGilli et al.l ||2007| ') model predicts 
a relatively large fraction of obscured AGN, ~50%, at high lu- 
minosities, while the observed value is only ^20%. The dashed 
gray line shows the e xpected dependen c e for a radiation-limited 
torus as described by IHonig fc Beckertl ||2007I V while the dotted 
gray line s hows the expectation for the original receding torus of 
ILawrenc3 119911) . both normalized to the observed value in the 



lO'* 



erg s 



bin. 



model is ~40% lower than the ITreister fc Urrvl (|2005l ) 
calculation, while it is now clear from recent Chan- 
dra and XMM observation s ()De Luca fc Molendil [20041 : 
iHickox fc Markevitchll2006D that a larger XRB intensity 
at low energies is more appropriate. 

In addition, in Figure [T2l we compare the observed de- 
pendence of the fraction of obscured sources on luminos- 
ity with the expectations for different geometrical param- 
eters of the obscuring material. If the height of the torus 
is roughly independent of luminosity, the change in cover- 
ing fraction is due to a change in inner radius (the orig- 
inal "receding torus" model), hence a rough L~^/^ de- 
pend e nce for the cont rast should be expected (jBarvainisI 
Il987t lLawrenc"i I1991D . If the effects of radiation pres- 
sure are incorpo r ated, in the case of a clumpy torus, 
iHonig fc Beckertl ()2007f ) derive a i"^/"* dependence for 
the contrast. As can be seen in Figure [T^l a de- 
pendence is too steep compared to observed data. This 
implies that the height of the obscuring material cannot 
be independent of the source luminosity and provides ev- 
idence for a radiation-limit e d stru cture, as the one sug- 
gested by IHonig fc Beckertl ()2007D . 

The dependence of the fraction of obscured AGN 
on redsh ift is more cont roversial. While some stud- 
ies (e.g., iLa Franca et al.|[2005; .Ballantyne et al.l 120061 : 
ITreister fc Urrvl 120061: iDella Ceca et al.l I2OO8D found a 



small increase in the fraction of obscured AGN at higher 
redshifts, o ther results sugge s t that this fraction is con- 
stant (e.g., [uidaitin [lOM lAkHiiliEin 11^^ The 
upper panel of Fig. |T3| shows the observed fraction of 
obscured AGN as a function of redshift for the sources 
in the ECDF-S. This fraction is high, ~90%, at z<l, 
while at higher redshifts it decreases to ^^30-40%. As 
for the luminosity dependence, in order to increase the 
significance of our results we added the E CDF-S sources 
to the sample of ' Treister fc Urrvl ()2006l ). The results 
for the large sample are consistent with those from the 
ECDF-S alone; namely, at low redshift there is a large 
fraction of obscured sources with a steep decline at Z'^l. 
This decline in the observed fraction of obscured sources 
can be easily explained by a simple observational fact: 
in order to be included in this sample a source needs 
to have a measured redshift from optical spectroscopy. 
For obscured AGN, t he optical light is dominated by 
the h ost galaxy (e.g., iBarger et al.ll2005l : ITreister et al.1 
120051 ). which becomes too faint for spectroscopy with 8 
meter class telescopes at Z'^l. In order to quantify this 
effect, we use the ratio of identified to total X-ray sources 
in a given optical magnitude bin (upper panel, fig. H]). 
This magnitude-dependent ratio is used to calculate the 
expected obscured AGN fraction including the X-ray and 
optical flux limits and the luminosity dependence of the 
obscured fra ction. This calcii l ation is done as described 
in detail by ITreister fc Urrvl (I2006D. Briefly, we use d 
the AGN population synthesis of lTreister fc Urrvl (|2005( l. 
toge ther with the library of AGN spectra described by 
Treis ter et al.l ([2004), and calculated the effective area of 
the survey as a function of both X-ray flux and optical 
magnitude, taking into account the spectroscopic incom- 
pleteness at each optical flux. This procedure outputs a 
expected number of observed obscured and unobscured 
AGN as a function of redshift, for a non-evolving ob- 
scured AGN fraction. 

As can be seen in the upper panel of Figure [T51 the ex- 
pected decrease is actually steeper than what is observed. 
This implies that the inferred fraction of obscured AGN 
in the ECDF-S should increase with redshift, once the 
selection effects in the sample are accoun ted for. This de- 
pende nce is consistent with the results of lTreister fc Urrvl 
(|2006l ). who found that this increase can be represented 
as a power-law of the form (l+z)", with a=0.4±0.1. 
These results are roughly independent of the method 
used to classify AGN; we obtained the same dependence 
using both our mixed X-ray/optical classification scheme 
and one based completely on optical spectroscopy. 

4.4. Number Density and Evolution of AGN 

Using our combined sample of X-ray selected AGN 
we can also compute the AGN spatial density as a 
function of redshift, and compare with the expecta- 
tions from existing AGN luminosity functions based 
on smaller samples. In order to calculate the AGN 
spatial density, we started from our collected sam- 
ple of 1377 sources described above. We then sep- 
arated the sources in low {Lx — ^0'^^^^'~*^ erg s^^). 



medium {Lx—W 



— in43-44.5 



erg s ) and high luminosity 



{Lx—W*^'^~'^^ erg s^"'^) bins. For each luminosity class 
we binned the sample in redshift so that each bin has at 
least 50 sources (20 for the highest luminosity sources). 
Then, the spatial density of each bin was calculated by 
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Fig. 13. — Fraction of obscured AGN as a function of red- 
shift. Upper panel: direct measurements using the sources on the 
ECDF-S field on ly (black circles with dashed error bars), from the 
ITreister fc Urrvl |^0()) sample (gray circles) and combining both 
samples (black circles with solid error bars). The expected ob- 
served fraction for an intrinsic fraction of 3:1 obscured to unob- 
scured AGN including optical and X-ray selection effects is shown 
by the black solid line. As can be seen, while the observed fraction 
of obscured AGN declines toward higher redshifts, if the X-ray and 
optical selection effects and the luminosity dependence of the ob- 
scured AGN fraction are taken into account this decline should be 
even stronger. Bottom panel: Inferred fraction of obscured AGN 
relative to an intrinsically constant fraction after correcting for 
selection effect and including the luminosity dependence of the ob- 
scured AGN fraction. Symbols are the same as for the upper paneh 
The corrected fraction of obscured AGN increases with redshift fol- 
lowing a power-law of the form (l+z)" with q=0 .4iIi0.1, consistent 
with the results found by ITreister fc Urrvl i200^ . 

summing the values of for each source in the bin, 
where Vc is the total comoving volume per unit area in 
that bin multiplied by the area covered by our super- 
sample at the X-ray flux of the source. In addition, upper 
limits for the spatial density on each bin were calculated 
by multiplying the values of Vc by the fraction of spec- 
troscopically identified sourc es at the optical magnit ude 
of the AGN, as described in ITreister fc Urrvl (|2006f l for 
the general sample and in the upper panel of Fig. [2] for 
the ECDF-S. 

In Fig. [M] we show the observed AGN spatial den- 
sity as a function of redshift for low, medium and high 
luminosity sources compared with the expected density 
obtained from integra ting the hard X-ray luminosity 
functi ons pr esented by lUeda et a l. (2003), Barge i~e"t al.l 
(|2005f l and iLa Franca erall (pOOSf ). While the X-ray 
sources in these studies were selected in similar ways, 
th e samples have dif ferent n umbers of sources ( 247 AGN 
inlUeda et ;ir'2003. 746 in iBarger et"aI1 120051 . and 508 
in iLa Franca ct al. 200^ with the emphasi s set at dif- 
feren t flux levels (mostl y bright sources in |Ueda et al.l 
l2003l mo derate fluxes in iLa Franca et aT]|2005l and faint 
sources in lBarger et al.ll2005l ). In addition, the modehng 
of the luminosity function is also different in these papers. 
While Bargcr ct al. (20dll) assumed pur e luminosity evo- 
lution , both Ucda ct al. (2003) and La Franca et al.l 
()2005[ ) found better results using a luminosity-dependent 
density evolution. As can be seen in Fig. [141 at high and 



Fig. 14. — Space density of AGN as a function of redshift de- 
rived from the sources in the ECDF-S. The filled circles show the 
observed density versus z in three luminosity bins, log Lx=41.5-43, 
43-44.5 and 44.5-48 erg s~^. Arrows show the effects of correcting 
these values for optical incompleteness effects, while the empty tri- 
angles show the results using the N'^^^/N^'^^ method to correct 
for incompleteness in the X-ray data. The Unes show the AGN 
space density derived using the luminosity functions of Ueda et 
al. (2003; solid lines), Barger et al. (2005]^ dashed lines) and La 
Franca et ah (2005; dotted lines). The |Barger ct al. (2005:) lu- 
minosity function is only defined for 2;<1.2. While in the highest 
luminosity bin all the studied luminosity functions agree very well 
with the observations, at intermediate luminosities a good agree- 
ment is found only at z<l. In the lowest luminosity bin, the effects 
of incompleteness are more important, and hence only after cor- 
recting for them a good agreement with the lum i nosity function of 
Ueda et al. ( 2003) is found. The IBarger et aP l(200a 'l luminosity 
function was not corrected for incompleteness. 

moderate luminosities all the luminosity functions stud- 
ied here agree well with the observations at low redshifts, 
while at high red shift s the expectations from both the 
lUeda et al.l (|2003D and lLa Franca et all (|2005D works are 
significantly above the observed values for intermediate 
luminosity sources. 

This discrepancy can most likely be explained by the 
effects of incompleteness in the X-ray data, which are 
particularly important for the lo wer luminosity sources. 
As it was concluded by e.g., iLa Franca fc Cristianil 
(|1997f ). the method used here is particularly sen- 
sitive to the effe cts of incompleteness. Th e work of 
lUeda et al.l (|2003D and ILa Franca et al.l (|2005D both used 

the A^o^s/iV™"^^ method, which attempts to account for 
incompleteness by using the expected number of sources 
in a given luminosity and redshift bin. One obvious 
caveat of this method is that the results are model de- 
pendent. We attempt to correct for incompleteness in 
the X-ray data, in the two lowest luminosity bins, where 
these effects are more important. In order to do that, 
we use the model described in the previous section in 
order to calculate the fraction of sources missed by the 
X-ray selection in a given luminosity and redshift bin due 
to the effects of obscuration. For the lowest luminosity 
sources, this correction is roughly ~2x, while at higher 
luminosities it is about 30-40%. We did not attempt to 
correct the highest luminosity bin, since the effects of 
absorption are negligible for these sources. As can be 
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seen in Fig. I14[ after the correction for incompleteness is 
applied a good agreement with the pr edicted spatial den - 
sity using the luminosity function of lUeda et al.l ()2003l ) 
is found. The remaining discrepancy with the work of 
iLa Franca et al.l ((2005|) could be explained because in 
this case even highly obscured Compton-thick AGN up 
to 7Vg=10^^ cm~ ^ are included, while in the case of 
lUeda et al.l (|2003D and in our work only mildly Compton- 
thick AGN with 7Vjy<10^^ cm~^ are considered. Finally, 
it is worth m entioning that in the luminosity function of 
iBarger et al.l ((2005) not attempts were made to correct 
for incompleteness in the X-ray selection. Hence, a good 
agreement is found only for the uncorrected data points. 

4.5. Infrared (IR) to X-ray ratio 

Emission at X-ray wavelengths, in particular in the 
hard band, is often used as a tracer for the direct AGN 
output mainly from Compton-scattered accretion disk 
photons, while the luminosity at longer wavelengths, in 
particular in the mid-far infrared, is associated with the 
re-radiation of the energy absorbed by the surrounding 
gas and dust (e.g., I Pier fc KrolikI IT9931 ) . There is cur- 
rently a strong debate about the geometry and charac- 
teristics of this surr ounding dust, in par ticular whether it 
has a smooth (e.g. .[Pier &: Kroliklll992f ) or clumpy (e.g., 
iKrolik fc Begelman|ll988D distribution. The ratio of IR 
to X-ray luminosity can be used to distinguish between 
thes e two distribution s and to constrain the dust geome- 
try (jLutz et al.ll2004D . For example, smooth torus mod- 
els predict large differences in the IR to X-ray ratio for 
obscured and unobscured sources, due to the significant 
effects of self-absorption. Contrarily, radiation transfer 
models of clumpy dust torii predict very small or no dif- 
ferences between the IR to X-ray ratio for obscured and 
unobscured AGN. 

The o bservational e videnc e remains controversial. The 
work of iHorst et al.l (|2006D reports that no differences 
were found in the mid-IR (measured at a fiducial rest- 
frame wavelength of 12.3 /zm) to X-ray ratio for a sample 
of 17 nearby AGN observed with the VLT-VISIR, which 
provides a relatively high angular resolution of ^^ ^ 0.35" . 
In apparent contradiction, iRamos Almeida et al.l ()2007f ) 
found a slightly smaller value for the X-ray to nuclear 
mid-IR (at 6.75 /xm) ratio from a sample of 57 AGN. 
For the latter, the observations were carried out using 
the ISOCAM camera onboard ISO, with a more limited 
angular resolution of ~4". iHorst et al.l (|20()8[ ) argue that 
precisely this difference in angular resolution explains 
the discrepant results. Using their high spatial resolu- 
tion images, they claim that the contribution from star 
formation, unresolved in the ISOCAM observations, can 
account for the observed differences between obscured 
and unobscured AGN. On the othe r hand, the relatively 
small sample of IHorst et al.l (|2006D traces larger intrin- 
sic luminosities for unobscured sources compared to the 
obscured sample, thus contributing to explain why no 
difference in the IR to X-ray ratio was found. Given the 
large scatter in the IR to X-ray ratio reported by both 
groups, a large sample of sources is required to reach 
statistically significant conclusions. 

In Figure [TSl we present the values of the ratio of L\ to 
Lx as a function of A, the rest-frame wavelength. These 
values were computed using both the IRAC and MIPS 
fluxes, for the sources in our sample with Spitzer detec- 



tions and Nh measurements. We also added the sources 
in the central re gion with Nh mea sured from the X- 
ray spectrum bv iTozzi et al.l ()2006D . In order to cor- 
rect the hard X-ray (2-8 keV) luminosities, Lx, for the 
effects of absorption we used t he photoel ectric absorp- 
tion c ross sections derived by iMorrison fc McCammonI 
(|1983D and the Nh values measured from X-ray spec- 
tral fitting as described above. Given that the effects 
of X-ray absorption in sources classified as unobscured 
are typically very small, and those values can be signifi- 
cantly affected by measurement errors, we only corrected 
the X-ray luminosities of the obscured sources. If no cor- 
rection for absorption is done to the X-ray luminosities 
of the obscured sources, our conclusions are unchanged. 
As can be seen in Fig [151 sources classified as obscured 
and unobscured AGN have significantly different aver- 
age values of L\/Lx, and this separation depends on 
wavelength. At the shortest wavelengths, A^l ^m, the 
separation in this ratio is rather small, which can be ex- 
plained by the contribution of the stellar light of the host 
galaxy to the integrated emission. This effect remains 
visible until A~2 ^m, where the stellar light starts to 
fade and emission from the host dust in the inner region 
of the surrounding material in the AGN begins to domi- 
nate. This emission is highly affected by self-absorption, 
as iriost torus models predict (see, e.g., Jr'ig r fc KroliH 
[1991 iNenkova et al.ll200a lH5nig etall l2006f) . thus ex- 
plaining why at these wavelengths unobscured sources 
have significantly higher values of Lx/Lx- At longer 
wavelengths, A>10 /im, the contrast between obscured 
and unobscured sources is reduced again, since the op- 
tical depth is reduced at longer wavelengths. Thus, it 
is expected that at A~30-40 /im the IR emission should 
become isotropic again. 

In order to compare with the observed averages for 
nearby Seyfert galaxies, in Figure [TSl we also present 
the composite IR spectrum for local sources classified 
as Seyfert 1 (unobscured) and S eyfert 2 (obscured), as 
compiled bv lPolletta et al.l (|2007f ). normalized at a fidu- 
cial wavelength of 100 /im, where the IR re-emission is 
expected to be fully isotropic. It is remarkable that both 
composite spectra agree well with our average values, 
thus indicating that the sources in our sample, most of 
them at z=0.5-l have very similar IR spectra to local 
active galaxies and higher luminosity sources. In ad- 
dition, th e results presented in Figure [151 allow to ex- 
jlain why IRamos Almeida et al.l p007f ) and IHorst et al.l 



(| 2008[ ) reach apparently discrepan t conclusions. While 
the work of IRamos Almeida et al.l (|2007f ) was based in 
observations at ~6 /xm, where the contrast between 
obscured and uno bscured sources is nearly maximal, 
IHorst et al.l ()2008[ ) used observations of a limited sam- 
ple of sources at ~12 /im, where the contrast is smaller 
than at short er wavelengths. T hus, it is possible that 
the results of IHorst et al.l (|2008f) are dominated by the 
intrinsic scatter in this relations, in particular given the 
low number of sources in their sample. 

In Figure [151 we also compare the observed infrared 
luminosities with t he predictions from t orus models re- 
cently presented bv INenkova et al.l (|2008D . These models 
assume a clumpy distribution for the obscuring material. 
While establishing the physical parameters of the dust 
surrounding the central region is a difficult task, which 
is beyond the scope of this paper, in Fig. [151 we show 
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Fig. 15. — Monochromatic luminosity normalized by the hard 
X-ray luminosity as a function of rest-frame wavelength for ob- 
scured (small black circles) and unobscured (small light gray cir- 
cles) AGN. For each source the four Spitzer IRAC and the MIPS- 
24 fim photometric data points were used. The large circles with 
error bars show the average values for each AGN type, in bins that 
include 40 sources. The dotted lines show the composit e spectra 
for ob scured and unobscured quasars, as compiled bv Poll etta et al.l 
l|2007l ) . while the dashed lines show exam ples of expec ted IR s pec- 
trum from the clumpy torus models of .Nenkova et aljl ('2008'), as 
described in the text. A clear separation between the average val- 
ues for obscured and unobscured sources is observed at wavelengths 
~5-15 fj,m, which can be explained by the effects of self-absorption 
by the surrounding material in obscured sources. 

examples of the expected IR spectrum for obscured and 
unobscured sources. The assum ed parameters are de- 
scribed by Nenk ova et al.l ()2008f) in the caption of their 
Figure 4, for a r~^ radial distribution of clouds. Fol- 
lowing the basic assumption of the original AGN uni- 
fication paradigm, the only difference between the ob- 
scured and unobscured model is the viewing angle. As 
can be seen, a decent agreement is found between the 
model spectrum and observations, in particular at wave- 
lengths longer than ~5 /xm. At shorter wavelengths, the 
additional contribution from the AGN host galaxy (not 
included in the model spectrum) starts to be significant, 
in particular for the obscured sources. 

Taking advantage of the large number of sources in 
our sample, in Figure fT6l we show L\/Lx as a function 
of wavelength for sources separated in three luminosity 
bins: Lx<10'*3 gj.g g-i^ ^O'^^ gj.g s-i<Lx<10''^ erg s'^ 
and Lx>10^* erg s~^. A clear conclusion from this fig- 
ure is that the difference in the average ratio between 
obscured and unobscured sources is largest at the lowest 
luminosity bin, while for higher luminosity sources the 
L\/ Lx ratios become very similar. In t heir study of the 
mid-IR properties of nearby AGN, Poll etta et~aLl ()2007f ) 
found somewhat consistent results. According to their 
Figure 9, the ratio ijj^ (at rest-frame 6 /im) to Lx for 
obscured sources, obtained combining their AGN2 and 
SF classes, is significantly lower than the ratio for unob- 
scured AGN (their AGNl class). In our interpretation of 
Figure 1161 this increasing contrast for lower luminosity 
sources can be understood in terms of the luminosity de- 
pendence of the geometrical parameters of the absorbing 



Fig. 16. — Same as Fig. llSl but separating the sample in three lu- 
minosity bins: logLx=42-43, 43-44 and >44 erg s~^. The largest 
separation between obscured and unobscured sources is found at 
the lowest luminosities, while for more luminous sources this dif- 
ference is smaller and almost negligible. One way to interpret this 
result is incl uding a depend ence of the opening angle on luminosity, 
as found by ITreister et al.l ||2008| ) , such that the lower luminosity 
sources have larger amounts of self-absorption due to a smaller 
opening angle. 

region, in par t icular the opening angle, as described by 
ITreister et all ()2008D . 

5. CONCLUSIONS 

We presented here the first results from our optical 
spectroscopy program, aimed to provide identifications 
and redshifts for the X-ray sources detected by Chan- 
dra in the ECDF-S field. As part of this work we tar- 
geted 339 X-ray sources from the catalog of IVirani et al.l 
(|2006f l and obtained redshifts and identifications for 160 
sources. We also added 26 sources with identifications 
obtained from the literature. While most of the sources 
with broad emission lines are found at z>l, the mean 
redshift for the sources in our sample is 0.75. In order 
to separate obscured and unobscured AGN we adopted a 
mixed scheme, combining the optical spectroscopy infor- 
mation with the indication of the X-ray spectrum given 
by the hardness ratio. This hybrid scheme agrees very 
well in most cases with the amount of neutral Hydrogen 
column density measured in the X-ray spectrum for the 
brightest sources. 

The optical colors of the obscured and/or low lumi- 
nosity AGN are dominated by the host galaxy, and are 
on average redder than non-active galaxies, as was found 
previously by other authors. We also focused on the 
fraction of obscured AGN and its possible dependence 
on luminosity and redshift. For the latter, we con firmed 

at higher significance the results of [Trcistcr fc Urrvl 

(|2006i) . indicating that after correcting for selection ef- 
fects the fraction of obscured AGN increases with red- 
shift. Regarding a possible luminosity dependence, we 
confirmed the functional form derived bv lTreister fc Urrvl 
(2005 ) and foun d that the population synthesis model of 
iGilli et al.l (j2007 ) significantly overestimates the fraction 
of obscured AGN at high luminosities. We also com- 
puted the AGN spatial density as a function of redshift 
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and found results in agreement wi th the expected va l- 
ues using the luminosity function of lBarger et alJ (|2005f ). 
but sig nificant discrepanc ies w ith the luminosity func - 
tions of lUeda et all (|2003[ ) and iLa Franca et all (|2005D . 
which can however by most likely explained by the effects 
of incompleteness in our X-ray-selected sample. 

Taking advantage of the deep multiwavelength data 
available in the ECDF-S we studied the infrared proper- 
ties of the X-ray selected AGN with spectroscopic iden- 
tifications. We found a significant difference in the frac- 
tion of bolometric light emitted at mid-IR wavelengths 
by obscured and unobscured AGN. These differences can 
be explained by the effects of dust self-absorption, for 
which the maximum contrast should be at ~5-15 ^m. 
Furthermore, by separating our sample into luminosity 
bins we found that the contrast is larger for the lower 
luminosity sources. One possible interpretation for this 
result is that the opening angle is larger for high lumi- 
nosity sources, such that the effects of self absorption are 
less important. 
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TABLE 1 

Log of Spectroscopic Observations 



Run 


Instrument 


Masks 


Slit Width 


Avg. Seeing 


Total Sources 


X-ray Sources 


Efficiency'' 


26-27/10/2003 


IMACS 


4 


1.0" 


~1" 


291 


74 


57% 


4-7/2/2005 


IMACS 


2 


1.2" 


~0.8" 


180 


64 


57% 


2-3/11/2005 


IMACS 


2 


1.2" 


~0.6" 


194 


29 


66% 


25-27/10/2006 


IMACS 


3 


1.0" 


0.5-1.5" 


280 


143 


32% 


21-22/11/2006 


IMACS 


1 


1.2" 


~0.8" 


109 


17 


20% 


18-20/2/2007'" 


IMACS 


1 


1.2" 


~1" 


109 


17 


20% 


Period 78 


VIMOS 


4 


1.0" 


1" 


283 


96 


54% 



^ Defined as the fraction of X-ray sources identified. The mask from the Nov. 2006 run was re-observed in order to 
obtain ~5 hours of integration time. 



TABLE 2 

X-RAY, OPTICAL AND MID-IR PROPERTIES OF THE X-RAY SOURCES IN THE ECDF-S 



ID=' 


X-ray Flux (> 


erg cm^s 


Hardness Ratio 


Redshift 


Instrument 


Class. <i 


log (X-ray Lum.) 


24 fim (ni 


>Jy) 




r 




Nh 




Softb 


Hard<= 


Value 


Upper 


Lower 








(erg s-i) 


Flux 


Error 


Source*^ 


Value 


Upper 


Lower 


(cm~^) 


2 


2.040e-15 


4.690e-15 


1.000 


0.000 


0.000 


9.990^ 


IMACS 


UNK 




51.27 


1.42 


F 










3 


2.670e-15 


6.620e-15 


-1.000 


0.000 


0.000 


9.990 


IMACS 


UNK 




223.50 


1.33 


F 


1.900 






20.843 


4 


3.650e-14 


6.590e-14 


-0.354 


0.030 


0.031 


2.011 


IMACS 


OAGN 


45.273 


229.90 


1.60 


F 


1.836 


0.123 


0.117 


20.941 


6 


6.760e-15 


8.590e-15 


-0.498 


0.069 


0.075 


3.031 


VIMOS 


UAGN 


44.821 


482.20 


1.22 


F 


2.045 


0.402 


0.347 


21.087 


7 


1.260e-15 


4.300e-15 


-0.054 


0.126 


0.146 


0.289 


IMACS 


ELG 


42.044 


92.09 


1.26 


F 











Note. — This table is published in it s entirety in t he ele ctronic edition of the Astrophysical Journal. A portion is shown here for guidance regarding its form and content. 

X-ray ID number from the catalog of IVirani et al.l 1120061) '' 0.5-2 keV"" 2-8 keV'^ UNK: Unkown, OAGN: Obscured AGN, UAGN: Unobuscred AGN, ELG: Emission line galaxy, ALG: Absorption lino 
galaxy, STAR: star.^ Survey in which the 24 //m source is identified. F: FIDEL, G: GOODS. ^ A 5:— 9.990 corresponds to a targeted source for which the redshift could not be measured. 
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